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Abstract 

We discuss experiments and observations aimed at testing the possible space-time variability of 
fundamental physical constants, predicted by the modern theory. Specifically, we consider two of 
the dimensionless physical parameters which are important for atomic and molecular physics: the 
fine-structure constant and the electron-to-proton mass ratio. We review the current status of such 
experiments and critically analyze recent claims of a detection of the variability of the fine-structure 
constant on the cosmological time scale. We stress that such a detection remains to be checked by 
future experiments and observations. The tightest of the firmly established upper limits read that 
the considered constants could not vary by more than 0.015% on the scale ~ 10^" years. 
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1 Introduction 

Possible variability of fundamental physical constants is one of the topical problems of contemporary 
physics. The modern theory (Supersymmetric Grand Unification Theory - SUSY GUT, Superstring 
models, etc.) has established that the coupling constant values which characterize different kinds of 
interactions (i) are "running" with the energy transfer and (ii) may be different in different regions of the 
Universe and vary in the course of cosmological evolution (e.g., Ref. [Q). The energy dependence of the 
coupling parameters has been reliably confirmed by high-energy experiments (see, e.g., Ref. [^), whereas 
the space-time variability of their low-energy limits so far escapes detection. 

Note that a numerical value of any dimensional physical parameter depends on arbitrary choice of 
physical units. In turn, there is no way to determine the units in a remote space-time region other than 
through the fundamental constants. Therefore it is meaningless to speak of a variation of a dimensional 
physical constant without specifying which of the other physical parameters are defined to be invariable. 
Usually, while speaking of variability of a dimensional physical parameter, one implies that all the other 
fundamental constants are fixed. So did Milne |^ and Dirac Q in their pioneering papers devoted 
to a possible change of the gravitational constant G. More recently, a number of authors considered 
cosmological theories with a time varying speed of light c (e.g., Ref. ||^). However, if we adopt the 
standard definition of meter |^ as the length of path traveled by light in vacuum in 1/299 792 458 s, 
then c = 2.997924 58 x 10^" cms~^ identically. Similarly, one cannot speak of variability of the electron 
mass me or charge e while using the Hartree units (ft = e = rrie = 1), most natural in atomic physics. 
Thus, only dimensionless combinations of the physical parameters are truly fundamental, and only such 
combinations will be considered hereafter. 

At present, the most promising candidate for the theory which is able to unify gravity with all other 
interactions is the Superstring theory, which treats gravity in a way consistent with quantum mechanics. 
All versions of the theory predict existence of the dilaton - a scalar partner to the tensorial graviton. Since 
the dilaton field (j) is generally not constant, the coupling constants and masses of elementary particles, 
being dependent on ^, should vary in space and time. Thus, the existence of a weakly coupled massless 
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dilaton entails small, but non-zero, observable consequences such as Jordan-Brans-Dicke-type deviations 
from General Relativity and cosmological variations of the gauge coupling constants ||^ . These variations 
depend on cosmological evolution of the dilaton field and may be non-monotonous as well as different in 
different space-time regions. 

In this paper, we focus on the space-time variability of the low-energy limits of two fundamental 
constants which are of paramount importance for atomic and molecular spectroscopy: 

(i) the fine-structure constant a = /he (Sommerfeld parameter), 

(ii) the electron-to-proton mass ratio /i = me/mp (Born-Oppenheimer parameter). 

The next section presents a compendium of the basic methods allowing one to obtain restrictions on 
possible variations of fundamental constants. In Sects. || and ||, we consider recent estimates of the values 
of a and ^, respectively, at cosmological redshifts z = 1-4 which correspond to epochs ~ 7-13 billion 
years ago. Conclusions are given in Sect. H 



2 Tests of variability of fundamental constants 

Techniques used to investigate time variation of the fundamental constants may be divided into extra- 
galactic and local methods. The latter ones include astronomical methods related to the Galaxy and the 
Solar system, geophysical methods, and laboratory measurements. 



2.1 Local tests 

2.1.1 Laboratory measurements 

Laboratory tests are based on comparison of different frequency standards, depending on different com- 
binations of the fundamental constants. Were these combinations changing differently, the frequency 
standards would eventually discord with each other. An interest in this possibility has been repeatedly 
excited since relative frequency drift was observed by several research groups using long term compar- 
isons of different frequency standards. For instance, a comparison of frequencies of He-Ne/CH4 lasers, 
H masers, and Hg+ clocks with a Cs standard |[ |l^, 11 1 has revealed relative drifts. Since the con- 



sidered frequency standards have a different dependence on a via relativistic contributions of order a^, 
the observed drift might be attributed to changing of the fine-structure constant. However, the more 
modern was the experiment, the smaller was the drift. Taking into account that the drift may be also 
related to some aging processes in experimental equipment, Prestage et al. |]lT| concluded that the current 
laboratory data provide only an upper limit \a/ce\ < 3.7 x 10^^"* yr~^. 



2.1.2 Analysis of the Oklo phenomenon 

The most stringent limits to variation of the fine-structure constant a and the coupling constant of the 
strong interaction have been originally inferred by Shlyakhter from results of an analysis of the 
isotope ratio ^''^Sm/^^^Sm in the ore body of the Oklo site in Gabon, West Africa. This ratio turned out 
to be considerably lower than the standard one, which is believed to have occurred due to operation of 
the natural uranium fission reactor about 2 x 10^ yr ago in those ores. One of the nuclear reactions in 
the fission chain was the resonance capture of neutrons by -'^^^Sm nuclei. Actually, the rate of the neutron 
capture reaction is sensitive to the energy of the relevant nuclear resonance level Er, which depends 
on the strong and electromagnetic interaction. Since the capture has been efficient 2 x 10^ yr ago, in 
means that the position of the resonance has not shifted by more than it width (very narrow) during 
the elapsed time. At variable a and invariable as (which is just a model assumption), the shift of the 
resonance level would be determined by changing the difference between the Coulomb energies of the 
ground-state nucleus ^''^Sm and the nucleus ^^°Sm* excited to the level Er. Unfortunately, there is no 
experimental data for the Coulomb energy of the excited ^^°Sm* in question. Using order-of-magnitude 
estimates, Shlyakhter ||l^ concluded that \a/a\ < 10~^^ yr^^- From an opposite model assumption that 
Us is changing whereas a ^constant, he derived a bound \as/cts\ ^ 10~^^ yr~^- 
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Damour and Dyson 1 13 performed a more careful analysis, which resulted in the upper bound \a/a\ < 
7 X 10~^^ yr"^- They have assumed that the Coulomb energy difference between the nuclear states of 
^^^Sm and ^^°Sm* in question is not less than that between the ground states of ^^^Sm and ^^°Sm. The 
latter energy difference has been estimated from isotope shifts and equals w 1 MeV. However, it looks 
unnatural that a weakly bound neutron (w 0.1 eV), captured by a ^'*^Sm nucleus to form the highly 
excited state ^^"Sm*, can so strongly affect the Coulomb energy. Moreover, heavy excited nuclei often 
have Coulomb energies smaller than those for their ground states (e.g., Ref. |l4|). This indicates the 
possibility of violation of the basic assumption involved in Ref. p^ , and therefore this method may 
possess a lower actual sensitivity. Furthermore, a correlation between a and as (which is likely in the 
frame of modern theory) might lead to considerable softening of the above-mentioned bound, as estimated 
by Sisterna and Vucetich |13] . 



2.1.3 Some other local tests 

Geophysical, geochemical, and palcontological data impose constraints on a possible changing of various 
combinations of fundamental constants over the past history of the Solar system, however most of these 
constraints are very indirect. A number of other methods are based on stellar and planetary models. The 
radii of the planets and stars and the reaction rates in them are influenced by values of the fundamental 
constants, which offers a possibility to check variability of the constants by studying, for example, lunar 
and Earth's secular accelerations. This was done using satellite data, tidal records, and ancient eclipses. 
Another possibility is offered by analyzing the data on binary pulsars and the luminosity of faint stars. 
Most of these have relatively low sensitivity. Their common weak point is the dependence on a model of 
a fairly complex phenomenon, involving many physical effects. 

An analysis of natural long-lived a- and /3-decayers in geological minerals and meteorites is much 
more sensitive. For instance, a strong bound, \a/a\ < 5 x 10~^^ yr""'^, was obtained by Dyson |l^ from 
an isotopic analysis of natural a- and /3-decay products in Earth's ores and meteorites. 

Having critically reviewed the wealth of the local tests, taking into account possible correlated syn- 
chronous changes of different physical constants, Sisterna and Vucetich [^sj derived restrictions on 
possible variation rates of individual physical constants for ages t less than a few billion years ago, 
which correspond to cosmological redshifts z < 0.2. In particular, they have arrived at the estimate 
a/a = (-1.3 ± 6.5) x IQ-^^ yr'^. 

All the local methods listed above give estimates for only a narrow space-time region around the Solar 
system. For example, the epoch of the Oklo reactor (1.8 x 10^ years ago) corresponds to the cosmological 
redshift z « 0.1. These tests cannot be extended to earlier evolutionary stages of the Universe, because 
the possible variation of the fundamental constants is, in general, unknown and may be oscillating [ pT| , |^. 
Another investigation is needed for higher cosmological redshifts. 



2.2 Extragalactic tests 

Extragalactic tests, in contrast to the local ones, concern values of the fundamental constants in distant 
areas of the early Universe. A test which relates to the earliest epoch is based on the standard model of 
the primordial nucleosynthesis. The amount of ^He produced in the Big Bang is mainly determined by 
the neutron-to-proton number ratio at the freezing-out of n<->p reactions. The freezing-out temperature 
Tf is determined by the competition between the expansion rate of the Universe and the /3-decay rate. 
A comparison of the observed primordial helium mass fraction, Yp — 0.24 ± 0.01, with a theoretical 
value allows one to obtain restrictions on the difference between the neutron and proton masses at the 
epoch of the nucleosynthesis and, through it, to estimate relative variation of the curvature radius R 
of extra dimensions in multidimensional Kaluza-Klein-like theories which in turn is related to the a 
value ijl^, [T^. However, as noted above, different coupling constants might change simultaneously. For 
example, increasing the constant of the weak interactions Gp would cause a weak freezing-out at a lower 
temperature, hence a decrease in the primordial '*He abundance. This process would compete with the 
one described above, therefore, it reduces sensitivity of the estimates. Finally, the restrictions would 



Comments on Atomic and Molecular Physics, 2000 



4 



be different for different cosmological models since the expansion rate of the Universe depends on the 
cosmological constant A. 

The most unambiguous estimation of the atomic and molecular constants at early epochs and in distant 
regions of the Universe can be performed using the extragalactic spectroscopy. Accurate measurements 
of the wavelengths in spectra of distant objects provide quantitative constraints on the variation rates of 
the physical constants. This opportunity has been first noted and used by Savedoff pO| , and in recent 
years exploited by many researchers (see, e.g., Refs. ||2^, and references therein). At present, the 
extragalactic spectroscopy enables one to probe the physical conditions in the Universe up to cosmological 
redshifts z < A, which correspond, by order of magnitude, to the scales ^ 10^*^ yr in time and ^ 10^ 
parsec in space. In the following sections, we review briefly the studies of the space-time variability of 
the fine-structure constant a and the electron-to-proton mass ratio /i, based on the latter method. 



3 Non-variability of a 

We have already mentioned in Sect. |2.1.l| that several laboratory tests hinted at a tentative time variation 
of a, but were later refuted by measurements at a higher level of accuracy. A similar situation has occurred 
for extragalactic tests at larger space-time scales. 

Bahcall and Schmidt were the first to use spectral observations of distant quasars to set a bound 
on the variability of the fine-structure constant. They have obtained an estimate Aa/a = (—2 ±5) x 10~^ 
ai z — 1.95. Later statistical analyses of fine-structure doublet lines in quasar spectra appeared to 

indicate a tentative variation of a (of the order of ^ 0.3% at the cosmological redshift z ~ 2). However, 
this tentative variation has been shown to result from a statistical bias . 

Another statistical examination of the fine-doublet wavelengths of absorption lines in quasar spectra 
| |2l[ | indicated a tentative (at the 2-3cr level) variability of a values by ~ 0.1% over the celestial sphere (as 
function of angle) at redshifts z ~ 2-3. However, this result has not been confirmed by a later analysis 
| |2^ , which was based on higher-quality spectra and yielded an order of magnitude higher precision. 

Quite recently, Webb et al. |^ have estimated a by comparing wavelengths of Fe ii and Mg ii fine- 
splitted spectral lines in extragalactic spectra and in the laboratory. Their result suggests a time- variation 
of a at the incredibly high accuracy level of ^ 10^'^%: the authors' estimate reads Aa/a = (—1.9 ± 
0.5) X 10~^ at z = 1.0-1.6. Note, however, two important sources of a possible systematic error which 
could mimic the effect: (a) Fe ii and Mg ii lines used are situated in different orders of the echelle-spectra, 
so relative shifts in calibration of the different orders can affect the result of comparison, and (b) if the 
isotopic composition varies during the evolution of the Universe, then the average doublet separations 
should vary due to the isotopic shifts. Were the relative abundances of Mg isotopes changing during the 
cosmological evolution, the Mg ii lines would be subject to an additional z-dependent shift relative to the 
Fe II lines, quite sufficient to simulate the variation of a (this shift can be easily estimated from recent 



laboratory measurements 1 28 ) 



The method based on the fine splitting of a line of the same ion species is not affected by these two 
uncertainty sources. We have studied the fine splitting of the doublet lines of Si iv, C iv, Mg ii and other 
ions, observed in spectra of distant quasars. According to quantum electrodynamics, the relative splitting 
of these lines SX/X is proportional to (neglecting small relativistic corrections, recently estimated by 



Dzuba et al. |g9| ) . We have selected the results of high-resolution observations |3ll |3^ , most suitable 
for an analysis of the variation of a. According to our analysis, presented elsewhere |?3[, the most reliable 
estimate of the possible deviation of the fine-structure constant at z — 2-4 from its present (z = 0) value: 

Aa/a = (-4.6 ± 4.3 [stat] ± 1.4 [syst]) x 10"^ (1) 

Thus, only an upper bound can be derived at present for the long-term variability of a: 

\a/a\ < 1.4 X 10"^'' yr'^ (2) 



(at the 95% confidence level). 
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4 Non-variability of /i 

The dimensionless Born-Oppenheimer constant /i = m,f./mp approximately equals the ratio of the con- 
stant of electromagnetic interaction a — jhc ~ 1/137 to the constant of strong interaction ols = 
jhc ~ 14, where g is the effective coupling constant calculated from the amplitude of 7r-meson-nucleon 
scattering at low energy. 

An early limit on the possible variation of this constant, Im/mI < 1-2 x 10~^° yr~^, has been derived 
from the concordance of K-Ar and Rb-Sr geochemical ages . The first astrophysical bound , based 
on the agreement between redshifts of atomic hydrogen and other lines in quasar absorption spectra, 
turned out to be twice stronger. 

Orders-of-magnitude more precise analysis has become possible due to discovery of a system of 
H2 absorption lines in the spectrum of quasar PKS 0528—250 at z = 2.811. A study of this system yields 
information about physical conditions and, in particular, the value of /x at this redshift (corresponding 
to the epoch when the Universe was several times younger than now). A possibility of distinguishing 
between the cosmological redshift of spectral wavelengths and shifts due to a variation of arises from the 
fact that the electronic, vibrational, and rotational energies of H2 each undergo a different dependence 
on the reduced mass of the molecule. Hence comparing ratios of wavelengths \i of various H2 electron- 
vibration-rotational lines in a quasar spectrum at some redshift z and in laboratory (at z = 0), we can 
trace variation of [i. We have calculated ||2l|, ^ sensitivity coefficients Ki of the wavelengths with 
respect to possible variation of /x and applied a linear regression analysis to the measured redshifts of 
individual lines Zi as function of Ki. If the proton mass in the epoch of line formation were different from 
the present value, the measured zi and Ki values would correlate: 

We have performed a z-to-K regression analysis using a modern high-resolution spectrum of PKS 
0528—250. Eighty-two of the H2 lines have been identified. The resulting parameter estimate and la 
uncertainty is 

A^Ji/^Ji = (-11.5 ± 7.6 [stat] ± 1.9 [syst]) x 10"^ (4) 
The 2a confidence bound on A/i//i reads 

|A/i//i| < 2.0 X lO^'^. (5) 

Assuming that the age of the Universe is ^ 1.5 x 10^° yr the redshift of the H2 absorption system 
z = 2.81080 corresponds to the elapsed time w 1.3 x 10^" yr (in the standard cosmological model). 
Therefore we arrive at the restriction 

l/i/Ml < 1.5 X lO"^'' yr"^ (6) 
on the variation rate of /it, averaged over 90% of the lifetime of the Universe. 



5 Conclusions 

Despite the theoretical prediction that fundamental constants of Nature should vary, no statistically 
significant variation of any of the constants has been reliably detected up to date, according to our 
point of view substantiated above. The upper limits obtained indicate that the constants of electroweak 
and strong interactions did not significantly change over the last 90% of the history of the Universe. 
The striking tightness of these limits is really astonishing and has already ruled out some theoretical 
models (see Refs. |2^, |2^, Q). A more elaborated theory (e.g., Ref. |7|) cannot be ruled out yet, but its 
parameters can be severely restricted (e.g., see Ref. [^2|). This shows that more precise measurements and 
observations and their accurate statistical analyses are required in order to detect the expected variations 
of the fundamental constants. 
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